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We study the effects of substructure on the rate of dark-matter annihilation in the Galactic 
halo. We use an analytic model for substructure that can extend numerical simulation results to 
scales too small to be resolved by the simulations. We first calibrate the analytic model to numerical 
simulations, and then determine the annihilation boost factor, for standard WIMP models as well as 
those with Sommerfeld (or other) enhancements, as a function of Galactocentric radius in the Milky 
Way. We provide an estimate of the dependence of the gamma-ray intensity of WIMP annihilation 
as a function of angular distance from the Galactic center. This methodology, coupled with future 
numerical simulation results can be a powerful tool that can be used to constrain WIMP properties 
using Fermi all-sky data. 
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I. INTRODUCTION 

Weakly-interacting massive particles (WIMPs) provide 
perhaps the most promising class of dark-matter candi- 
dates. These are particles that arise in theories of new 
electroweak-scale physics, such as low-energy supersym- 
metry (SUSY) U or models with universal extra dimen- 
sions (UEDs) [2]. Experiments that seek to directly de- 
tect these particles in low-background experiments, or to 
indirectly detect them through observation of energetic 
neutrinos from WIMP annihilation in the Sun, are now 
beginning to dig into the favored WIMP parameter space. 
However, there are also prospects for indirectly detect- 
ing WIMPs through observation of gamma rays and/or 
cosmic-ray positrons, antiprotons, or antideuterons from 
WIMP annihilation in the Galactic halo. These annihi- 
lation products have received considerable attention in 
the recent literature with the attribution of some re- 
ported (and still controversial) cosmic-ray anomalies [3] 
to WIMP-annihilation products. 

The total rate at which WIMPs annihilate in the 
Galactic halo is proportional to the volume integral of 
the square of the dark-matter density. In the canonical 
(and simplest) model, the dark matter is smoothly dis- 
tributed in the halo in a spherically-symmetric way with 
a dark-matter density p{r) that is a monotonically de- 
creasing function of r; for example, it is often modeled as 
an isothermal or a Navarro-Frenk- White (NFW) profile 

However, analytic arguments and numerical simula- 
tions show that there should be substructure in the dark- 
matter distribution in the Galactic halo [5-7 . The dark 
matter may be clumped; some of it may be bound in 
higher-density self-bound subhalos; and some may be in 
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tidal streams [5] . This is an outcome of hierarchical clus- 
tering, in which small dense halos form first and then 
merge to form progressively larger structures. If dark 
matter in the halo is clumped, then the total annihila- 
tion rate will be enhanced by some boost factor [S], the 
increased rate per unit volume in dense regions outweigh- 
ing the decreased rate per unit volume in lower-density 
regions. 

In the canonical WIMP scenario, the substructure may 
have a roughly scale-invariant distribution in substruc- 
ture mass/size extending all the way down to substruc- 
tures on mass scales ~ 10"^° Mq [71[Tni[II], roughly 22 
orders of magnitude smaller than the ~ 10"'^^ Mq Milky 
Way halo. While state-of-the-art numerical calculations 
now have the resolution to simulate several decades in 
this hierarchy, they are very far from being able to follow 
the survival and evolution of the smallest substructures 
through all of the generations in the structure-formation 
hierarchy that result in a Milky Way halo. Analytic 
calculations of these survival probabilities are difficult 

P,rT3]. 

However, given that the smallest subhalos are likely to 
be the densest, the boost factor may depend significantly 
on the existence of these substructures. This is partic- 
ularly true in models with a Sommerfeld enhancement 
in the annihilation rate |14j — those where the annihila- 
tion rate increases with lower WIMP velocities — since the 
smallest subhalos are also likely to have the smallest ve- 
locity dispersion. 

In a previous paper [15J, we presented an analytic 
model to describe the self-similar substructure expected 
from hierarchical clustering. The model predicted a high- 
density power-law tail for the probability for a given 
point in the halo to be in a clump of density p. Sub- 
sequent to that paper, another appeared [T5] presented 
N-body simulations that showed this power-law tail. Af- 
ter calibration to numerical simulations, the model can 
be used to extrapolate the results of numerical simula- 
tions to substructure-mass scales far smaller than those 
currently accessible to the simulations. We investigated 
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in Ref. the dependence of the boost factor assuming 
a canonical WIMP-annihilation rate (i.e., no Sommerfeld 
enhancement). A complementary analytical approach to 
this problem, based on the stable clustering hypothesis, 
was recently presented by [T7] . 

In this paper we re- visit and extend those calculations. 
We first describe in Section [IT] the analytic model. We 
extend the model by including a finite width for the 
smoothly distributed component of dark matter in ad- 
dition to the high-density power-law tail. We then use in 
Section |III| state-of-the-art N-body simulations to mea- 
sure the width of the smooth component and the am- 
plitude of the power-law tail and thus calibrate the an- 
alytic model to current simulations. We moreover de- 
termine how the substructure distribution varies with 
Galactocentric radius in the halo. In Section ITVl we use 
the calibrated model to determine the boost factors for 
canonical WIMPs and for WIMPs with a Sommerfeld 
enhancement to the annihilation rate. Section |V] re- 
views the substructure-model parameters expected from 
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WIMP models, and Section VI determines the angular 
dependence of the intensity of gamma-ray radiation both 
with and without the substructure boost factor we ob- 
tain. Section IVIII reviews the model and results and then 
provides some comments, caveats, and directions for fu- 
ture development of the model. 



II. SUBSTRUCTURE MODEL 

Before discussing our analytic model for substructure 
we begin by presenting the canonical smooth-halo model 
against which the substructure model will be compared. 
We take as the canonical model for the Galactic halo an 
NFW profile, 
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as a function of Galactocentric radius r, with parameters 
Ps = 0.051 GeV cm~^ and = 21.7 kpc taken to provide 
a reasonable fit to the Milky Way rotation curve. 

If dark matter in the halo is clumped, then the densi- 
ties at all points with the same r will not necessarily 
be the same. Instead, there will be some probability 
distribution function P{p,r) defined so that P(p,r)dp 
is the probability that a particular point in the Galac- 
tic halo (at some fixed Galactocentric radius r) has a 
density between p and p + dp. According to the argu- 
ments of Ref. '_15j, a fraction fs of the volume of the 
halo (in fact most of the volume; fs — 1, as we will see) 
should be filled with a smooth dark-matter component 
with density ph , and a fraction 1 — <^ 1 will consist of 
a high-density clumped component, with something like 
a power-law distribution of densities. Thus, the proba- 
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where 8(x) is the Heaviside step function. Here, the first 
term describes the smooth host halo component as having 
a Gaussian distribution in In p with mean density ph and 
a Gaussian width in Inp of A. This distribution peaks 
at a density po — e^'^'^ ^^Ph, a density slightly below p/j. 

The second term in Eq. ^ is the high-density power- 
law tail due to substructures that remain from earlier 
generations in the structure-formation hierarchy. The 
parameter ph in Eq. ^ is in fact a function of r; we will 
see below that Phir) ~ p to very good accuracy. There 
may also be further r dependence in P{p; r) through an 
r dependence of the parameters A, a, and fa- As we will 
see, our simulations are not yet good enough to allow the 
r dependence of a to be resolved, and so we take it to 
be a constant. The simulations are, however, sufficiently 
resolved to see a strong dependence of /^(r) on r, which 
we detail below. 

The distribution P(p; r) in Eq. ^ is normalized so 
that P{p] r) dp — 1. It can be integrated to give the 
mean density. 
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This thus provides a relation between the mean density 
Phi^) of the smooth component and the canonical host 
halo density p. Strictly speaking, this relation requires 
the maximum density Pmax out to which the power-law 
tail extends. However, we will see that numerically 1 — fs 
is so small that Phif) — p{r) is a very good approximation 
for practical purposes. 



III. CALIBRATION TO SIMULATIONS 

We make use of the Milky Way-like halo in the Via 
Lactea II simulation [18] and calibrate our model to the 
P(p; r) measured from the simulated particle distribu- 
tion. Fig.[l]shows the PDF P{p) derived from Via Lactea 
II for three Galactocentric radii, at 10 kpc (lowest solid 
line, red), at 100 kpc (middle line, blue), and at 300 kpc 
(highest line, cyan). To obtain this P{p), we consider 
an ellipsoidal shell that follows an iso-density contour in 
the halo. For each particle in that shell, we calculate 
di = p/p where p is estimated in the usual manner from 
the nearest N neighbors using a symmetric SPH kernel 
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fraction is well approximated by 
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FIG. 1. The probability distribution function P{5) obtained 
from simulations. The solid curves are the simulation results 
at r = 300, 100 and 10 kpc (from top to bottom). The dashed 
curves show our analytic approximations to the power law tail. 
The dotted curve indicates the contribution to the finite width 
of the smooth component at 10 kpc from Poisson fluctuations 
due to the use of A'^ = 32 neighbors in the density estimator. 
Note that the mean host halo density ph to which the x axis 
is normalized is ~ 175x (~ 4900 x) smaller for the 100-kpc 
(300-kpc) curve than in the 10-kpc curve. 



and the median density p is obtained from all the parti- 
cles in that ellipsoidal shell. These 5i are then binned in 
equally-spaced bins in log]^o((5). In each of these bins, we 
calculate Pilogi^S) — X^j'^i"^ where the sum is over all 
particles in that bin; the 6^^ weighting gives a volume- 
fraction distribution. The distribution in logj^g ^ then 
converted to a distribution in 5 and normalized. 
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Note that at radii less than ^ 20 kpc, 1 — /s(r) drops 
faster than Eq. CI); for example, 1 — /^(lOkpc) = 4 x 
10-4 « 1.5xl0-^p(10kpc)/p(100kpc))-°-26. This close 
to the center, however, the dumpiness of the simulated 
halo is likely artifically suppressed due to finite resolu- 
tion effects. The best-fit values of a are 0.0 ± 0.1 at all 
radii greater than 20 kpc. In the following, we implicitly 
assume a = and the radial dependence in /s given by 
Eq. 0. 
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FIG. 2. The probability distribution function P{5) at 
100 kpc for particle densities estimated from the nearest 
N = (16, 32, 64, 128, 1024) neighbors. 



A. Power-Law Tail 

The central features of Fig. [T] relevant here are the 
high-density power-law tails predicted by Ref. [TS] (and 
seen already in simulations [TB]). The figure shows that 
the amplitude of the high-density power-law tail is larger 
at larger radii. This can be attributed largely to the fact 
that the mean density p is ^ 175 times lower at 100 kpc 
than at 10 kpc, and another factor ~ 30 times lower at 
300 kpc, and so the ratio of the density in substructures 
to the mean density is higher at larger radii. 

We now use this simulation to calibrate the analytic 
model at a variety of radii r, from 4 to 300 kpc. At each 
radius we fit for the power law parameters a and /«. We 
find that at radii greater than ^ 20 kpc, the smooth-halo 



B. Finite Width of the Smooth Component 

The simulation results shown in Fig. [l] show a finite 
width A for the smooth component. However, care must 
be taken as Poisson fluctuations due to the finite num- 
ber N of nearest neighbors in the density estimator will 
also contribute to the width. In Fig. [2] we show P{5) 
at 100 kpc for densities determined with N = 16, 32, 
64, 128, and 1024. The dotted curves indicated the ex- 
pected contribution to the width from Poisson fluctua- 
tions (and note that the true and Poisson widths should 
add in quadrature), which we obtained by running the 
density estimator on a randomly distributed sample of 
10^ particles. As N is increased, the width of the smooth 
component decreases, but not quite as fast as the Poisson 
fluctuations, and by A'^ = 1024 it is clear that the true 
width has been resolved to be about A ~ 0.2. At 10 kpc 
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(not shown here) it remains unresolved, and we conclude 
only that A < 0.2 at radii less than 100 kpc. 



IV. ANNIHILATION BOOST FACTOR 

We now consider the boost of the annihilation rate 
in a halo with substructure relative to the rate in the 
canonical smooth-halo model. 

The annihilation rate (per unit volume) at any point 
in the Galactic halo is 



r = (av) 
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where {av) is the annihilation cross section (times rel- 
ative velocity v, averaged over the velocity distribution 
of the halo, and is the WIMP mass; i.e. {av) — 
dv f (v) av , where f (v) (normalized to f{v)dv = 
1) is the WIMP pairwise velocity distribution at that 
point in the Galactic halo. For the canonical WIMP, av is 
approximately velocity independent at galactic-halo ve- 
locities, and so (av) = (crw), a constant, liav depends on 
velocity, then (av) may vary from one point in the halo 
to another due to the possible variation of f{v) with po- 
sition in the halo. 



A. Standard velocity dependence 

Let us first consider the canonical WIMP where the 
cross-section factor (av) is velocity-independent. In this 
case, the change in the velocity dispersion that may ac- 
company clumping is irrelevant. The total annihilation 
rate in some volume V is then proportional to a volume 
integral of the density squared — i.e., oc / p'^dV. The 
annihilation rate for a fixed total mass of dark matter 
enclosed within a given volume will thus be enhanced by 
a "boost factor" if the matter is not uniformly distributed 
within that volume. For the canonical WIMP, this boost 
factor B{r) is the ratio 
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where the second equality follows since P{p, r) ~ 
{l/V){dV/dp) (with V the volume in the halo). The 
quantity Pmax is the maximum density. We take that to 
be ^ 80 GeV cm~^, which is about five times the mean 
virial density of the rare first collapsed structures in the 
Universe (178 times the mean density of the Universe at 
z ~ 40); we will say more later about Pmax- 

In the standard WIMP model, where ((Ti;)=constant, 
the boost factor B{r) is given by (noting that p ~ ph) 



B{r) = fse^ 
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There is r dependence in this boost factor via Eq. Q 
in /s, via Eq. ([T]) in w p, and in principle also in A 
although it turns out to be negligible. 

There are two contributions to this boost factor: The 
first comes from the finite width A of the smooth 
component — i.e., Bs = Js& — and it depends very 
strongly on A. Given that we find A < 0.2 in our sim- 
ulations, we infer that the boost factor due to the finite 
width of the smooth component is no more than a few 
percent. 

We now turn to the second term, that due to substruc- 
ture, the central focus of this work. Given that Pmax ^ P 
in the Milky Way halo, the boost factor would be essen- 
tially independent of Pmax if a were a > 1. However, we 
find in the simulations that a < 1 in a Milky- Way-like 
halo, in which case the boost-factor results for the canon- 
ical WIMP do indeed depend on pmax- The fact that the 
integral in Eq. ([6) is dominated by the high-density end 
implies that the boost factor due to clumping is deter- 
mined primarily by the early-collapsing highest-density 
(and lowest-mass) substructures. 




FIG. 3. The density probability distribution function P{p) for 
an isolated NFW halo with concentration c = 10, truncated at 
rvir. The P{p) ~ behavior continues beyond pvir (vertical 
dashed line) to ~ p{rs) = 1/12 c^/ f{c) pvir- 



What is an appropriate value for Pmax? A first guess 
might be the virial density of the earliest-collapsing ha- 
los, Pvir = 178 pcrit(^c) ^ 16 GcV cm'^ (zj40)3 (for 
a collapse redshift Zc ^ 0). However, depending on 
its density profile, most of a halo's volume might have 
considerably higher densities. In Fig. |3] we show the 
density probability distribution P(p) = 1/VdV/dp for 
an isolated NFW halo. At low densities, in the out- 
skirts of the halo where p ^ r^^, the density proba- 
bility falls as p"^, matching the power law tail of the 
local density probability function (Eq. In the in- 
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nermost regions of the halo, where p ~ , we have 
P{p) ~ p~'^ ■ The transition between these two regimes 
occurs at Ps = 1/12 c^/f{c) p^-„, where c — Rvi-c/ts is the 
concentration of the halo and /(c) = ln(l + c) — c/(l + c). 
We set /Omax = Ps = 1/12 c^/f{c} Pvir, which then de- 
pends only on the concentration with which the earliest- 
collapsing halos are born. Numerical simulations [7 in- 
dicate low natal concentrations of c w 2 — 5, correspond- 
ing to /Omax ~ 1.5—11 Pvir- For dcfiniteness, we pick 
an intermediate value of c — 3.5 and Zc = 40, giving 
Pmax = 80 GeV cm-3. 

Close to the center, the clumped fraction 1 — /s, as 
determined from the numerical simulation, is so small 
that the boost factor remains close to unity. As the mean 
halo density p^ decreases with radius while the clumped 
fraction increases (cf. Eq. Q), the local boost factor 
grows considerably in the outer regions of the halo. Note, 
however, that the total luminosity of the halo does not 
increase in proportion to this local boost factor. The 
overall luminosity is dominated by radii < Ts, and the 
total boost from substructure within a radius R must be 
evaluated numerically,. 



B{<R) 



CB{r)p{r) 



(8) 



We show in Fig. |4] the differential and cumulative lumi- 
nosity boost factor as a function of radius for the Via 
Lactea II host halo, assuming a = and 1 — fs{r) as 
given m Eq. Q. The boost factor remains close to unity 
in the center and only reaches 1.5 at the Sun's distance 
of 8 kpiQ implying that if the WIMP annihilation cross 
section has the canonical (av) cx constant dependence 
on the velocity, the local boost from substructure is un- 
likely to provide the missing factor of 100-1000 needed 
to explain the cosmic-ray anomalies (see also p2]). The 
total luminosity of the halo, however, can be appreciably 
boosted by substructure. The cumulative boost B{< r) 
increases to ~ 17 at the virial radius. 

Before proceeding further, we note that if a — (and 
Pmax ^ p), then the boost factor is B{r) = -f (1 — 
fs){Pmax/p), an expression that is easily understood. The 
first term is simply the usual annihilation rate due to the 
smoothly-distributed dark matter. The second is that 
due to clumping. If the integrand in Eq. ^ is dominated 
by the high-density end, it implies that most of the anni- 
hilation in the clumped component is taking place in the 
smallest and densest subhalos. If so, then the annihila- 
tion rate, per unit volume, from substructure should be 
proportional simply to the spatial density of the subhalos 
(i.e., how many are there per unit volume), which itself 
should be proportional to the ratio of subhalo to host 
halo density, (pmax/p); this is consistent with Eq. ^ if 
a = 0. 
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^ It would be only 1.14 using the value of 1 
kpc, instead of Eq.ffl. 
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FIG. 4. The local substructure boost B{r) [solid) and the 
cumulative luminosity boost B{<r) (dotted), as a function of 
radius. 



B. Sommerfeld enhancement 

Suppose now that the annihilation cross section is such 
that the thermally-averaged cross section is av cx v~^, 
where v is the rms relative velocity for annihilating 
WIMPs. There will then be an additional enhance- 
ment in the annihilation rate since lower-mass subhalos 
will have smaller velocity dispersions. We account for 
this additional effect as follows: We first recall that the 
characteristic density of a first-generation halo collaps- 
ing at z = 40 with a concentration of c = 3.5 is = 
80 GeV cm^'^, and that the corresponding characteristic 
velocity is Vg = \JG M{< rs)/rs ~ 1.0 x 10^"^ km sec~^. 
We then note that a typical Milky- Way host halo of 
mass 2 x IO^^Mq and concentration c = 15 at z = 
has a corresponding characteristic density and velocity 
of Ps = 0.076 GeV cm'^ and Vs = 200 km sec^^ This 
thus suggests a rough scaling v cx p^^ ''^. We emphasize 
that this scaling is only meant to very roughly capture the 
relation between density and relative velocity of the DM 
particles. In reality there likely is no such simple one- 
to-one relationship between these two quantities, since 
regions with similar densities can be bound to subhalos 
of very different masses and hence have very different ve- 
locity dispersions. However, as will become clear below, 
our results are not very sensitive to the exact value of the 
power law exponent, as long as regions of higher density 
(at a fixed galacto-centric distance) typically have lower 
velocity dispersions. 

In the following, we consider only the Sommerfeld- 
enhanced boost factor from the clumped component and 
disregard the small contribution from the finite width 
of the smooth component. We further assume that the 
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dark- matter velocity dispersion wmw — 220 km sec~^ in 
the Galactic halo does not vary with Galactocentric ra- 
dius r. Strictly speaking, this constancy does not hold 
for a self-gravitating NFW distribution. Realistically, 
though, the Milky Way disk contributes very significantly 
to the potential in the inner Galaxy, and so the dark- 
matter spatial/ velocity distribution in the inner Galaxy 
cannot be a pure self-gravitating NFW distribution. Our 
assumption of a constant dark-matter velocity dispersion 
is probably closer to the truth than the radial change in 
the velocity dispersion implied by an NFW distribution. 

With these approximations and assumptions, the boost 
factor is 
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which we then integrate with the relation v[p) = 

'yo(/3/Pmax)"^ '^^ tO find 
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Again, for the values of a we see in the simulation, the in- 
tegrand in Eq. (|9| is dominated by the high end, increas- 
ingly so for /3 > 0. Now, even though 1 — /« may be small 
(the clumped fraction is small) , the velocity enhancement 
(umw/'^^o)^ may be large, even for values /3 ~ 1. Again, 
if the integrand is dominated by the high-density tail, 
it implies that most of the annihilation in the clumped 
fraction is occurring in the lowest-mass highest-density 
regions. And if so, then the annihilation rate per unit 
volume should again be proportional simply to the ratio 
of subhalo to host halo density (pmnx/p), which is again 
implied in Eq. (10) if a = 0. 



If the mass of the force carrier particle mediating 
the annihilation is non-zero, the Sommerfeld effect satu- 
rates at a finite velocity Usat, when the de-Broglie wave- 
length of the particle becomes longer than the range 
of interaction. With such a saturation the integral in 
Eq.Q is split into two parts: one from p to psat — 
PmJvo/vsatY^^''''^ with a velocity-dependent Sommer- 
feld enhancement term {vmw/'^{p))^ i ^-nd a second part 
from psat to Pmax with a constant enhancement factor of 

(-WMW/Wsat)'^, 

- + -/.)(! + «)f^)' 




Psat 
P 



-1.75/3^ 



if) 



(11) 



Fig. [5] shows the boost factor at the solar radius 
as a function of the velocity parameter /3. The solid 
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FIG. 5. The boost factor at the solar radius as a function 
of the parameter /3, for no saturation (solid), Vsat ~ 10~* c 
(dotted), and Usat ~ 10~^ c (dashed). The Sommerfeld-like 
enhancement grows as v''^ until Usat- 



line depicts the case without saturation, and the dot- 
ted and dashed lines with Wsat/c = 10^^ and 10^^, re- 
spectively. We use Pmax — 80 GeV cm~'^ and vq = 
1.0 X 10^'^ km sec~^ here. With Sommerfeld enhance- 
ment it is possible to get very large substructure boost 
factors even at the solar radius. This substructure boost, 
of course, applies in addition to the Sommerfeld enhance- 
ment of the smooth halo annihilation luminosity. 

Fig. [6] shows the cumulative boost factor B{<r) as a 
function of radius. In the top panel we plot curves for 
/3 = 0, 0.1, 0.3, 0.5, 0.7, and 1.0, assuming no satura- 
tion. In the bottom panel we fix /3 = 1 (except for the 
reference /? = case) and vary the saturation velocity, 
v^^t/c = 10"'*, 10-^ 10-6, 10"^, and 0. Increasing /3, or 
lowering Ugat at a fixed /3, leads to significant increases 
in the cumulative boost factor. For example, the total 
boost factor from within the virial radius of ^ 300 kpc 
grows from 17 (/3 = 0) to 1.3 x 10*^ for (3^1 without 
saturation. Even with Wsat/c — 10"^, the /? = 1 case still 
results in about an order of magnitude increase in the 
total boost, to ^ 120. 



V. APPLICATION TO WIMP MODELS 

We now assemble estimates for the numerical val- 
ues of Wo and Pmax for WIMP models. After freeze- 
out of WIMP annihilation in the early Universe, the 
WIMPs may continue to scatter from the more abundant 
light standard model particles. These scatterings sup- 
press perturbations in the WIMP density on sub-horizon 
scales until these elastic-scattering interactions cease; i.e.. 
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FIG. 6. The cumulative boost factor as a function of distance 
from the Galactic center. In both panels the solid curve rep- 
resents the canonical case (/3 — 0). Top panel: B{< r) for 
different values of /3 with no velocity saturation: /3 = 0.1, 0.3, 
0.5, 0.7, and 1.0 from bottom to top. Bottom panel: B{<r) 
for P = 1.0 and different values of the saturation velocity: 
Usat/c = lO""*, 10"^ 10"®, 10"'', and from bottom to top. 



at kinetic decoupling. This post-frcezeout kinetic cou- 
pling of WIMPs suppresses primordial perturbations on 
mass scales smaller than ~ 33(rkd/10 MeV)"^ Me 
[6], where Tkd is the kinetic-decoupling temperature. 
The smallest substructures in the Milky Way halo will 
therefore have formation masses no smaller than Ale- 
Close to the halo center, tidal interactions and im- 
pulsive stellar encounters may remove some of frac- 
tion of this mass [12;, but the dense cuspy cores con- 
tributing the majority of the annihilation luminosity 



likely survive [T3]. Detailed calculations of the rele- 
vant elastic-scattering reactions shows that in SUSY and 
UED models for WIMPs, this mass scale spans the range 
10"^ Me < Mc < 100 Me [U, the precise value de- 
pending on the particle-physics details. Objects in this 
mass range undergo gravitational collapse at a redshift 
Zc — 40 — logjQ (Mc/Me) HO]; the weak dependence of 
the collapse redshift on Mc arises from the n — >■ —3 
limit of the primordial power spectrum at small scales. 
These first collapsed objects obtain virial velocities t^vir = 
vo ~ 1.0 X 10-3 (Mc/Me)^/3(zc/40)^/^ km and virial 
densities pvir — 16 {zc/4:0)^ GeV cm-^, corresponding to 
Pmax = 80(ze/40)3 {c/3.5 ff{3.5)/f{c) GcV cm-3. 

For a given WIMP model, the cutoff mass Mc can be 
calculated and thus Pmax a-nd vq obtained. Given these 



parameters, our Eq. ( 10 1 can provide the boost factor, as 



a function of radius, for given Sommerfeld parameters (3 
and Usat- 

For example, suppose we would like a substructure 
boost of B{r) ~ 10'^ locally to account for reported 
cosmic-ray anomalies. For an Earth mass cutoff (Ale ~ 
Me, Vo = 1.0 X 10-3 km s-\ p^ax 80 GeV cm^^)^ ^^]^_ 
ing a = 1, and a Sommerfeld model without saturation. 



the boost factor in Eq. ( 10 ) becomes 



B{r) 



jl-fs) P, 
1 + 1.75/? 



nax / WmwV 



(12) 



which is > 10^ for /3 > 0.69. Likewise, assuming (3 — 1, 
a saturation velocity fsat/c < 3.3 x 10-^ is necessary to 
get a local substructure boost factor greater than 10"^. 



VI. ANGULAR DEPENDENCE OF THE 
GAMMA-RAY INTENSITY 



We now consider the dependence of the intensity (pho- 
tons cm-^ s-"'^ sr--'^) of gamma-ray radiation from WIMP 
annihilation in the Milky Way halo as a function of the 
angular separation -0 between a given line of sight and 
the Galactic center. This intensity /(i/') can be written 
as an integral 



dl 



P+Rl- 2lRo cos 



xB dP + Rt - 2lRocosij 



(13) 



along a line-of-sight distance I. If the halo is smooth, then 
B{r) — 1, and the integral is / oc J p"^ dl, of the square of 
the smooth-halo density along the line of sight. If all of 
the annihilation in the halo took place in highly dense and 
very small subclumps, then the intensity would depend 
on an integral J oc J pdl simply of the density (rather 
than density squared). If the intensity depends on the 
integral of p^, then the intensity will vary more rapidly 
with ip, rising rapidly toward the Galactic center, than if 
it depends on the integral of p, as shown in Fig. [7j 
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FIG. 7. Intensity of gamma-ray radiation from WIMP an- 
nihilation in the Milky Way as a function of the (cosine of 
the) angle that the line of sight makes with the Galactic 
center. All three curves are normalized to have the same in- 
tensity at cosip = 0. The top (solid) curve is the intensity 
due to annihilation in a smooth halo; i.e., for 7 cx J dl. 
The bottom (dotted) curve is that for annihilation entirely in 
subhalos; i.e., / oc J pdl. The middle (dashed) curve is for a 
Sommerfeld-enhanced annihilation [i.e., Eq. ( |11[ |] with /3 = 1, 
u,at/c = 5.0 X 10-^ 



Most generally, there may be annihilation in both the 
smooth component and in the clumped component, in 
which case the angular dependence of the intensity will 
fall somewhere in between , as illustrated in Fig. [t] 
Measurement of this angular variation may help shed em- 
pirical light on the existence of a boost factor. 



VII. DISCUSSION 

It has long been recognized that there may be a hierar- 
chy of structure in the Milky Way halo, with substructure 
existing all the way down to the mass scale of a fraction 
of an Earth mass, more than 22 orders of magnitude from 
the ^ 10^^ Mq Milky Way halo mass. It has also been 
noted that this substructure may have serious implica- 
tions for the dark-matter-annihilation rate in the halo. If 
the WIMP has the canonical velocity- independent (crw), 
then the boost factor could be as high as (pmax/p) ~ 200, 
if all of the substructure was preserved. If the WIMP has 
a Sommerfeld enhancement to the annihilation rate, the 
boost factor could be even higher, and perhaps dramati- 
cally so. 

Unfortunately, the 22 or more orders of magnitude be- 
tween the substructure cutoff mass and the Milky Way 
mass prevent the survival of the smallest-scale substruc- 
ture to be addressed directly with simulations, and reli- 



able analytic calculations of the survival fraction are sim- 
ilarly difficult. In earlier work |15j . we used the nearly 
self-similar behavior of hierarchical clustering to develop 
an analytic approach to estimate the substructure sur- 
vival. The central result of that work was a prediction 
that the probability P{p) for a given point in the Galac- 
tic halo to have a local density p will have a high-density 
power-law tail. 

Here, we have fit the parameters of that analytic model 
to new N-body simulations that can resolve the high- 
density tail in P{p). The analytic model then allows us 
to extrapolate the behavior of the simulations to mass 
scales far below the simulation's resolution scale. As a 
result, we have a simple analytic expression for the dis- 
tribution of dark-matter densities within the Milky Way 
halo, as a function of Galactocentric radius. The key 
qualitative result is that the fraction of the Milky Way 
volume occupied by substructures is small (1 — /s ^ 10~^ 
in the central regions, ^ 10^^ in the outskirts). In par- 
ticular, when we calculate the boost factor for WIMPs 
with the canonical (i.e., no) dependence of (av) on the 
velocity, we find that it is small: only about 50% at the 
solar radius, and only ~ 17 for the total boost within 
the virial radius. The PDF, supplemented with a scal- 
ing for the subhalo velocity dispersion with subhalo den- 
sity, allows us to also analytically estimate the boost 
factor from substructure in the presence of a Sommer- 
feld enhancement. The central result here is given in 



Eqs. (10 11), which provide a boost factor in terms of 
the power-law index /3 for the velocity scaling of (crw), 
the saturation velocity Wgat, the maximum substructure 
density Pmax — 80 GeV cm~'^, and the velocity dispersion 
Wo of the smallest halos. We estimate numerically that 
local (8 kpc) boost factors > 10^ can be obtained with 
/3 > 0.69, or for /3 = 1 with i;,at/c < 3.3 x 10"^. 

We also discussed the finite width of the density dis- 
tribution smoothly-distributed component of halo dark 
matter, but find that the effects of this finite width on 
the boost factor are small. 

There are several caveats to our results and several 
improvements that can be pursued in future work. (1) 
Hierarchical clustering is not a precisely self-similar pro- 
cess. In particular, given that the primordial mass power 
spectrum P{k) — >■ fc" has a power-law index that be- 
comes smaller at higher k (smaller distance/mass scales, 
earlier collapse times) , the power-law index a in our work 
will probably have some scale dependence, becoming, if 
anything, larger at higher densities. Taking this into ac- 
count, our estimates for the boost factors due to sub- 
structure are probably on the high side. (2) The small- 
scale mass cutoff Mc for canonical WIMPs TT may not 
apply if the new interactions required for a Sommerfeld 
enhancement are taken into account. The cutoff mass 
may therefore be very different. If it is much larger, than 
the boost factors will be reduced (substructures will not 
extend to such small scales). If the cutoff mass is much 
smaller, the boost factor may be increased relative to 
our estimates, but, given the weak dependence of Zc and 
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thus Pmax on Mc, not by much. (3) While identification 
of the high-density power-law tail in P{p) and measure- 
ment of its small amplitude is a big step forward, our 
N-body measurement of the parameters required to de- 
scribe the distribution P{p, r) can certainly be improved 
upon. It will be important in future work to measure 
fsir) more precisely, and to measure and determine the 
r dependence of a and A. (4) Fig. [s] indicates that at 
densities above Pmax, the a — power law for P{p) will 
steepen to something closer to a = 2. This steep high- 
density tail is due to the small-r p ^ 1/r dependence in 
the earliest NFW halos. For a canonical WIMP with a 
velocity-independent (av) , this steepening will not affect 
the results, since the integrand in Eq. ^ is proportional 
to and hence dominated by the low density (i.e., 
near /Omax) end of this steep tail. With a Sommerfeld 
enhancement, however, the integrand is proportional to 
pi.75 i3-a (jgj)-)^ which for a = 2 is dominated by the 
high-density end if /? > 8/7. In other words, if /3 > 8/7, 
the substructure annihilation enhancement will be dom- 
inated by the 1/r cusps in the earliest NFW subhalos. 



rather than the a = part of P{p) that we have consid- 
ered until now. A more detailed calculation of this effect 
would depend sensitively on the smallest radius at which 
the 1 /r NFW behavior is valid, and we leave such a study 
to future work. 
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